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■ Abstract: The dynamics of density and metric perturbations is investigated for the previously 

■ developed model where the decay of the vacuum energy into matter (or vice versa) is due to the 

O 



renormalization group (RG) running of the cosmological constant (CC) term. The evolution of 
the CC depends on the single parameter which characterizes the running of the CC produced 
^P' by the quantum effects of matter fields of the unknown high energy theory below the Planck scale. 

The sign of v indicates whether bosons or fermions dominate in the running. The spectrum of 
^ ' perturbations is computed assuming an adiabatic regime and an isotropic stress tensor. Moreover, 

' the perturbations of the CC term are generated from the simplest covariant form suggested by 

the RG model under consideration. The corresponding numerical analysis shows that for > 
there is a depletion of the matter power spectrum at low scales (large wave numbers) as compared 
to the standard ACDM model, whereas for v < Q there is an excess of power at low scales. We 
find that the LSS data rule out the range \v\ > 10~^ while the values \v\ < 10^^ look perfectly 
acceptable. For < the excess of power at low scales grows rapidly and the bound is more severe. 
From the particle physics viewpoint, the values ~ 10^^ correspond to the "desert" in the mass 
spectrum above the GUT scale Mx 10^^ GeV . Our results are consistent with those obtained in 
other dynamical models admitting an interaction between dark matter and dark energy. We find 
that the matter power spectrum analysis is a highly efficient method to discover a possible scale 
dependence of the vacuum energy. 
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1. Introduction 



The analysis of cosmic perturbations |||, ^, ^ represents one of the main tests for the candidate 
cosmological models and especially for the models of a time-dependent Dark Energy (DE). The 
recent data on the accelerated expansion of the Universe from the high redshift Type la supernovae 
Q, combined with the measurements of the CMB anisotropics |5|, ^, ^ and LSS ^ indicate that 
the DE is responsible for most of the energy balance in the present Universe. The natural candidate 
to be the DE is the cosmological constant, mainly if we take into account that from the quantum 
field theory (QFT) point of view the zero value of the CC would be extremely unnatural |jlO|, |l^ . 
The solution of the standard CC problems, such as the origin of the fine-tuning between vacuum and 
induced counterparts and the problem of coincidence, is not known. Historically the first attempts 



to solve some of these problems went along the lines of dynamical adjustments mechanisms |13] 



Nowadays the discussion is approached from many different perspectives and has generated an 
extensive literature, see e.g. |14, 15| and references therein. One of the most interesting aspects 



of these discussions is the possibility to have a slowly varying DE density of the Universe |16]. 



This option cannot be ruled out by the analysis of the present data [17, In case that the 
new generation of cosmic observational experiments will indeed detect a variable DE, the natural 
question to ask is: would this fact really mean the manifestation of a qualitatively new physical 



reality such as quintessence , Chaplygin gas |Q , extra dimensions ||2^ or of low-energy quantum 



gravity ||2^, 23, 24, 25 1? Let us emphasize here that one cannot honestly address this question 
without also exploring the possibility of a time- variable CC whose running is due to known physical 
effects, that is, due to the quantum effects of matter fields on a curved classical background p^ . 

The most natural manner to investigate the possible running of the CC produced by the 
quantum effects of matter fields is the renormalization group (RG) in curved space-time |^, |l^ . 
Unfortunately, the existing calculational methods are not sufficient for performing a complete 
theoretical investigation of this problem [^]. At the same time, the phenomenological proposal 
1 11, ^ shows that a positive result can not be ruled out. The hypothesis of a standard 



quadratic form of decoupling for the quantum effects of matter fields at low energies [11, 31] 



leads to consistent cosmological models of a running CC with potentially observable consequences 
1 32, 33, 34, |3^, 36 1. Among the many possibilities one can distinguish essentially two different 



models of running CC in QFT in curved space-time ^. The first one p^, 3C, admits that 



there is energy exchange between matter and vacuum sectors of the theory^. This leads to the 
vacuum energy density function /)a = Pk{H) associated to the cosmological term. A, while the 
Newton constant G remains invariable^. The second model does not admit a matter- vacuum 
energy exchange, but has both parameters scale-dependent G = G{H) and p\ = Pa{H) [^J. In 
the present work we shall concentrate on the analysis of the perturbations for the first model 
only, and postpone the investigation of the second model for another publication. Let us remark 



^Generalizations of this kind of semiclassical RG cosmological models (which may e.g. help to alleviate the cosmic 
coincidence problem) are also possible, see 

^See js^ for a general discussion of energy exchange between cosmic fluids. 

^Let us notice that here we follow |ll| and identify the cosmic energy scale with the Hubble parameter H, despite 



other choices are also possible [0 |l2||2|, |l|, | 



that the analysis of cosmic perturbations in the first model has already been addressed in |4C] 
using an indirect procedure based on bounding the amplification of the density matter spectrum 
in the recombination era caused by vacuum decay into CDM. This amplification is necessary to 
compensate for the dilution of the 5pm/ Pm matter spectrum at low redshifts resulting from the 
enhanced matter-radiation density pm ^ associated to the decay of the vacuum into CDM. A 
qualitatively similar analysis has been performed in ||4l[| , and also in |^2| using the so-called CMB 
shift parameter [^j. It is worth noticing that the upper bounds for the running obtained in all these 
references disagree between themselves, in some cases at the level of several orders of magnitude. 
This severe divergence illustrates the importance and necessity of performing a direct calculation 
of the perturbations, which we undertake here thoroughly. For example, the estimate used in 



|4C] was derived from the phenomenological study presented in Ref. [^], where the matter power 
spectrum of density fluctuations measured by the 2dF Galaxy Redshift Survey (2dFGRS) was 
compared with the corresponding matter power spectrum derived from the measurements of the 
CMB anisotropics. The comparison served to place a bound at the level of 10% on the maximum 
possible difference between the two derivations of the matter power spectrum 1 43 1. By requiring 
that the amplification of the density fluctuations of the RG model at the recombination era do not 
surpass this bound the authors of |4C] were able to obtain a corresponding bound on the the RG 
parameter of ||2^, specifically v < 10^^. Of course this bound is only approximate as it does not 
come from a direct calculation of the density perturbations in the RG model. Notwithstanding, 
as we shall see in what follows, there is a qualitative agreement between these results and the 
ones obtained in the direct calculation that we present here. In both cases the RG models with 
Planck-scale mass particles are not favored and one is led to conclude that a significant energy 
gap (or "desert") below the Planck scale, Mp, is necessary to avoid contradictions between the 
calculated spectrum of the density perturbations and the observational data. The favored mass 
scale turns out to be three orders of magnitude below Mp, and can be identified with the typical 
GUT scale Mx ~ 10^^ GeV . In this paper we will provide a rigorous quantitative derivation of 
this result. Let us also remark that the analysis of perturbations presented here is not covered by 
the many works in the literature dealing with perturbations in a Universe with dynamical dark 



energy - see e.g. [15] and references therein. In fact, the nature of the DE under consideration 



is not directly related to the properties of dynamical scalar field(s), but to the QFT running of 
the cosmological constant. Therefore, our approach is qualitatively new and has not been dealt 
with in the literature before. It constitutes, to our knowledge, the first computation of density 
perturbations in the presence of a running cosmological term. This is important in order to assess 
the impact on structure formation from a possible scale dependence of the vacuum energy. 

The paper is organized as follows. In the next section we summarize the necessary information 



about the variable CC model under discussion. The reader may consult the papers |11, 32, 33] 
for further technical details and extensive discussions. In section 3 we derive the equations for 
the perturbations. Section 4 is devoted to the numerical solution of these equations and to their 
comparison with the 2dFGRS H. Finally, in the last section we draw our conclusions. 



2. The background cosmological solution in the running CC model 



In the framework under consideration [0, 29|, with renormahzation group corrected CC and en- 
ergy exchange between vacuum and matter sectors, the cosmological evolution is governed by the 
following three ingredients: (i) Friedmann equation, 

H\z) = ^ [p^{z) + Pa{z)] + Hinl (1 + zf , (2.1) 

where Pm{z) and pk{z) are the densities of matter and vacuum energies taken as functions of 
the cosmological redshift, z, with HqVL^{1 + z)"^ = —kjc?\ (ii) Conservation law (following from 
the Bianchi identity) which regulates the exchange of energy between matter and vacuum, 

dpA ^ dprn ^ "iPm 2) 

dz dz 1 + 2 ' 

and {Hi) the renormahzation group equation for the CC density, assuming a soft form of decoupling 



for the cosmological constant |11, 2£, 31 1 



4 



"^^^ - ^ aH^M\ (2.3) 



d\nH (47r) 

In this equation, the effective mass M represents an additive sum of the contributions of all virtual 
massive particles, and cj = ±1 depending on whether fermions or bosons dominate at the highest 
energies. It proves useful introducing the new (dimensionless) parameter 

V = ^ . (2.4 

The value of v defines the strength of the quantum effects and has direct physical interpretation. 
For example, taking the "canonical" value ly = i^o = l/127r = 0.026 means that the effective 
cumulative mass M of all virtual massive particles is equal to the Planck mass, whereas |z^| ~ 1 
means the existence of a particle with trans-planckian mass or of nearly Planck-mass particles with 
huge multiplicities. Taking [z^l ~ 10~^ means that the spectrum of particles is bounded from above 
at the (GUT) scale Mx ^ W^^ GeV ~ IQ-^Mp. Much smaller values of < lO"'^ mean the 
existence of an extra unbroken symmetry between bosons and fermions (e.g. supersymmetry) at 
the GUT scale etc. 

We already noticed that Eq. (^]^) implies the possibility of an energy transfer between vacuum 
and matter, which may be interpreted as a decay of the vacuum into matter or vice versa. While 
the decay of the vacuum energy into ordinary matter and radiation could be problematic for the 
observed CMB, the decay into CDM can be allowed [16, ^ provided the rate is sufficiently small. 



It is of course part of the present study to evaluate the possible size of this rate. 



*For a dimensionless parameter this would be 0{p^ /m?) (that is the standard Appelquist and Carazzone 
decoupling law. In the special case of vacuum energy, which is a dimension-4 parameter, this leads to the quadratic 
dependence on the scale, which we associate to the Hubble parameter . Inclusion of bulk viscosity effects would 
lead to linear terms in H [UsI. We neglect them in the present work. 



The complete analytic solution of the equations (2.1), (2.2), (|2.3|) can be found in [p9|, 32] 



Here we present only those formulas which are relevant for the perturbations calculus. First of all, 
from (|2.3| ) it follows that 

p^{z-.)=pl + ^^Ml{H^-Hl). (2.5) 

Furthermore, we introduce notations for the pressure and density of the vacuum and matter com- 
ponents, such that the total values become (we assume pressureless matter at present) 

Pt = Pm + PK, Pt = Pm + Pa, where Pa = -PA , Pm = . (2.6) 
The quantities (^^) determine the dynamical evolution law for the scale factor: 

H + H^ = - = {pt + 3Pt) = —{pm- 2pA) . (2.7) 

ad 6 



What we shall need in what follows are the two ratios 



Pm _ {i + z){H^y -2Hini(i + 



i2 



2 



Pt 3[H^-Hlnl{l + z) 

, _ _ PA_ sH^-{i + z)W-Hlnl{i + zr 

where the prime indicates d/dz. The expansion rate H[z) in this model is given as follows: 

\ l — SiyJ I — u 1 — 6V 

In the last formula = J^5)M + svlth of the Dark Matter (DM) and baryonic matter 

densities (relative to the critical density) in the present-day Universe. Let us remark that the scale 
dependence of the vacuum energy (^.5|) does not imply a change of the equation of state for the 



"cosmological constant fluid". In this work we follow |29, 3C, 32 1 and assume that pA = —Pki 
even when quantum corrections are taken into account. It should be clear that this equation of 
state is valid irrespective of whether A is constant or variable. A different issue is that the models 
with variable cosmological parameters (like the present one) can be interpreted in terms of an 
effective equation of state corresponding to a self-conserved dark energy, as it is usually assumed 
in e.g. scalar field models of the DE. The general framework for this effective EOS analysis has 



been developed in |35, where a particular discussion of this issue for the present running CC 



model is also included. 



3. Deriving the perturbations equations 

In order to derive the equations for the density and metric perturbations we follow the standard 
formalism (see e.g. [Q]-[^]). We will perform our calculation of density perturbations under the 
following two assumptions: 



The entropic perturbations are negligible; 



• The full energy-momentum tensor is free of anisotropic stresses. 

In the absence of more information we believe that these hypotheses are reasonable as they aim 
at the maximum simplicity in the presentation of our results. Notice that in a single fluid case, 
the pressureless fluid and the vacuum equation of state have zero anisotropic pressure ; if both 
fluids interact, as it happens in the present case, a contribution can appear, but we cannot model 
it without introducing further free parameters. Similarly, the contribution of entropic perturba- 
tions must be introduced phenomenologically, and this would lead to the appearance of more free 
parameters, thus masking any constraint we can obtain on the fundamental parameter v of our 
RG framework. It should also be stressed that the same approximations are made in similar works 
in the literature, for example in reference where the density perturbations in the alternative 
framework of interacting quintessence models are considered and a corresponding bound on the 
parameter that couples dark energy and matter is derived. Therefore, in order to better compare 
our results with these (and other future) studies, and also to essentially free our presentation of 
unnecessary complications extrinsic to the physics of the RG running, in what follows we will 
adopt the aforementioned set of canonical hypotheses. At the same time we compute the density 
perturbations by neglecting the spatial curvature, which is irrelevant at the epoch when these 
linear perturbations were generated, as explained e.g. in . Nevertheless we have kept the spatial 
curvature in the background metric where it could play a role at the present time. 

Let us first introduce the 4-vector velocity [7^. In the co-moving coordinates C/'^ = [/q = 1 
and W = Ui = 0. The total energy-momentum tensor of matter and vacuum can be expressed via 
dp) as 

T; = [pt + Pt)U-U^-Pt5-^, (3.1) 

such that = pt and If = -Pt dj. It is pretty easy to derive the covariant derivative of [7^: 

V^U'' = 3H. (3.2) 

The last equation is important, for it enables us to perturb the Hubble parameter and eventually 
the running CC. In particular, using ( |3.2D we can rewrite Eq. ( |2.5| ) in the form 

PA = A + B{V^U'')\ (3.3) 

where A = pi - Mj, Hi , B = . 

In writting the cosmological term as in equation (^j^), we intend to obtain a covariant expres- 
sion which reproduces the background relations. Although this expression may not be unique 
(in the sense that additional terms vanishing in the background could nevertheless contribute to 
the perturbations), any other covariant form reproducing that behavior would be more involved. 
Therefore, we will adopt the simplest possibility ( |3.3| ) as our ansatz for the calculation. Any vari- 
able cosmological term that varies as the square of the Hubble function could be written in the 
same manner. However, the concrete way such term appears in the background, and the relation 
of the factors A and B to the quantum field mass through the parameter u ~ see Eq. (^.4|)- is 
specific of the present model. 



Consider next simultaneous perturbations of the two densities and the metric 



Pm ^ Pm (1 + <Jm) , PA ^ PA (1 + ^h) i 9liu ^ g'^j.u = 9fiu + . (3.4) 

The background metric g^j^y = diagjl, —a?{t)5ij^ corresponds to the solution described in the 
prrevious section. We assume the synchronous coordinate condition /iq^ = and obtain the 
variation of the i?Q in the form 

^^o° = -|^ftV^#f^V (3.5) 
It proves useful to introduce the notation 



dt \ a2 

This variable satisfies the equation 

h + 2Hh = SvrG {5pm - 25pa) ■ (3.7) 

Perturbing the covariant conservation law, V ^Ty = 0, we arrive at the following equation for the 
f = Q component: 

5pm + (6* - ■^) Pm + Si? 5pm = -6pA, (3.8) 

and, for the = i components, at the equation 

PmO + bH PraO + PmO = 5- (^PA • (3.9) 

Here we have used the constraint 5U^ = 0, and also the notation \Ii{5U'^) = 6 for the covariant 
derivative of the perturbed 3-velocity. We should not confuse the wave number k in Eq. ( |3.9D and 
below with the spatial curvature parameter introduced before, by the context it should be obvious. 
Indeed, it is understood that we have written all the previous perturbation equations in Fourier 
space, namely using the standard Fourier representation for all quantities: 

[■ f]3i. 

/(x, t) = —- /(k, t) e^^-^ , k = \k\. (3.10) 



J (2vr)3 

Perturbing the Eq. (|3.3|) and using the relation (3.2) one finds 



Pa-Sa = 6pA = 2b(9- ^) {V^Un = 6HB - |) • (3.11) 

It proves useful to rewrite the perturbations in terms of another set of variables. Let us 
introduce the new variable v = 0pm/pt = ^ /i and trade the time derivatives for the ones in 
the redshift parameter z, using the formula d/dt = —(1 + z)Hd/dz. After a straightforward 
calculation the perturbation equations can be cast into the following form: 

/ , (3/1-5) kHl + z)h , 
v+^^v = 5a, (3.12) 



/( 3/2 

/i 



1 



(l + z)// \/i 2 



(/2 5a)' 



3/2 
1 + 2 



(5a 



(3.13) 



h! 



2h 

1 + z 



2z/ 



{1 + z)q 



ifl (^m - 2 /2 (Ja) 



(3.14) 




Sa = 

In the last two equations we have introduced the function 

2iyH(z) 



(3.15) 



3HHz) - 3H^Ql{l + zY 



(3.16) 



The equation ( 3.15 ) is not dynamical, rather it represents a constraint which can be replaced in 
the other equations (3.12)-(3.14). Using notations ( p.g| ) and ( |3.16| ) we eventually arrive at: 



, (3/1 - 5) 

V H — V 

1 + z 



Hfi 



k'^Q{l + z) 



2H 



(3.17) 



h' + 



2{v - 1) _ 2u (2v _ h 
l + z (1 + z) Wi 



(3.18) 



5L + 



/( 3/2 
/i 



l + z 



1 



+ 



l + z 



1 



(3.19) 



These equations constitute the complete set of coupled Fourier modes for the velocity, density 
and metric perturbations in the presence of a running cosmological term. One can check that for 
vanishing u they reproduce the situation in the ACDM model |]l|, ||, |^, as expected from the fact 
that = corresponds to having a strictly constant pA and conserved matter density pm- The 
previous set of equations are in the final form to be used in our numerical analysis of the density 
perturbations for arbitrary values of v and their impact on the large scale structure formation in 
the present Universe^. 



4. Power-spectrum analysis and comparison with the galaxy redshift survey 
4.1 Power spectrum and transfer function 

In the previous section we have obtained the coupled set of differential equations ( |3.17 ), ( p. 18 ) 
and ( |3.19| ) describing the dynamics of the density and metric perturbations. In order to perform 
the numerical analysis of these equations, we have to fix the initial conditions and also define 
the limits for the variation of the redshift parameter z and for the perturbations wave number 



^In Ref. perturbations were considered in the very early Universe for an inflationary model with exponentially 
decaying vacuum energy under the assumption of a thermal spectrum and with no obvious relation to QFT and RG. 



k = |k|. The analysis performed here appUes after the radiation dominated era, in fact z must 
vary from about z = 1100 (the recombination era) to z = (today). For the sake of the numerical 
analysis we measure k in the units of /iMpc~^, where h is the reduced Hubble constant. In 
these units let us notice that we must consider k < 0.15 /iMpc^^ since the observational data 
concerning the linear regime are in this range . Expressing the scales in terms of the Hubble 
radius Hq^ =~ 3000 /iMpc implies that we must consider the values of k between and A50Hq. 

The expression for the baryon matter power spectrum is given in many places of the litera- 
ture [Q, ^. We use the concrete form given in [^]. At z=0 we have^ 

where = Clpj + 17^, with = rJ^JM + defined before. (Notice that for flat geometry 

ily = 1.) Here A is a normalization coefficient (see below), T[k) is the transfer function and g{^) 
is the growth function ~ see below for details on these functions. We consider a scale invariant 
(Harrison-Zeldovich) primordial scenario as revealed by the assumed linear dependence in /c of (4.1). 



The transfer function takes into account the physical processes occurring around the recombination 
era. The use of the transfer function allows to perform the numerical analysis and the comparison 
with the recent data from the 2dFGRS |Q deep in the matter dominated era, for example in the 
range of z between zq = 500 and Zf = 0. The choice of zq is quite arbitrary, the only requirement 
is that the system could evaluate sufficiently deep into the matter dominated phase. Essentially 
the results do not depend on the precise value of zq. 

In order to fix properly the initial conditions, we use the transfer function presented in refer- 
ences [51, 4^. This transfer function assumes a scale invariant primordial spectrum, and determines 



the spectrum today considering the Universe composed of dark matter and a cosmological constant 
that generates the actual accelerated expansion phase. With the aid of this transfer function we 
can fix the initial conditions at a redshift after the recombination. Of course, our model contains 
the traditional cosmological constant as a particular case when our relevant RG parameter z/ = 0, 
see Eq. (U). 



Let us now specify the transfer function. To this end we use the following notations [49|: 

r = n%he-^'h-{^%/^lr), g = qik)= ' (4.2) 

[hT) Mpc 



In these expressions, T is Sugiyama's shape parameter j50|] . Given an initial spectrum (e.g., 
Harrison-Zeldovich), the power spectrum today can be obtained by integrating the coupled system 
of equations for the evolution of the Universe and the Boltzmann equations. The transfer func- 
tion for the baryonic power spectrum can be approximated by the following numerical fit (BBKS 
transfer function) 



^(^) = ^"olT!f^^ + 3-89g + (16.1g)' + (5.46g)=^ + (6.71g)^l . (4.3) 



''In general the power spectrum depends on a parameter n labeling the possibility of a non-vacuum initial state 
of the perturbations of quantum origin We take of course n = 0, corresponding to the power spectrum of the 



Furthermore, the growth function appearing in (|4.1|) reads |52| 



(4.4) 



This function takes into account the effect of the cosmological constant in the structure formation (a 
suppressing effect when the CC is positive) . Notice that both the transfer function and the growth 
function are involved in the computation of the matter power spectrum in Eq. (^]^). Finally, the 
normalization coefficient A in (4.1) can be fixed using the COBE measurements of the spectrum of 
CMB anisotropies. This coefficient is connected to the quadrupole momentum Qrms 

of the CMB 



anisotropy spectrum by the relation |4£] 



A 



(2/^)4 ^ ^ 

^ ^ 5 r2 



(4.5) 



where Ih = ^ SOOO/i"^ Mpc is the present Hubble radius and Tq = 2.725 ± 0.001 is the 
present CMB temperature. The quadrupole anisotropy will be taken as 



Qr 



(4.6) 



This value is obtained from COBE normalization and is consistent with the more recent results 



of the WMAP measurements |53] using the prior of a scale invariant spectrum. Taking all these 
estimates into account, we can fix 



yl = 6.8 X 10^ Mpc"^ , 



(4.7) 



which corresponds to the initial vacuum state for the density perturbations used in |49]. 

There is evidently nothing sacred about this particular choice of the transfer function. Many 
others have been proposed in the literature. For example, the following useful simplified transfer 
function has been proposed by Peebles in Ref. Q: 

4.7 P 



T{k) 



1 + 



8k 



+ 



(4.8) 



In this function it is understood (as in the BBKS one (|4.3| )) that k is given in /iMpc~^ units (as 
mentioned above). Strictly speaking the normalization used in ( [4.^ ) should be re-computed for 
the present model using the CMB anisotropy spectrum, but this calculation is very involved and 
lies beyond the scope of the present work. Let us notice that the remaining arbitrariness in the 
normalization of the power spectrum may not affect the qualitative results of our investigation. In 
particular, we have checked that the results obtained with Peeble's function ( [4. 8] ) are very close 
to those obtained with the BBKS function ( |4.3| ). However, the last one is more detailed, and we 
can control better the contribution of each component of the matter content for the final power 



spectrum. For this reason we present here only the results for the transfer function (4.3) 



4.2 Initial conditions. Normalization with respect to to the ACDM model 

Let us now explain in some more detail how do we fix the initial conditions for our specific problem. 
The strategy is the following. We first consider the standard ACDM cosmological constant case 



as a reference to normalize our calculation. In this way the baryonic spectrum today (that is, for 
z = 0) can be described by e.g. the transfer function ( |4.3| ). This allows to fix the conditions for 
the density contrast 5m Sit z = for the ACDM case. Moreover, for the initial condition on the 
metric function, h, we suppose that it is identical to the density contrast, because for the ACDM 
model this is so (up to a small factor). Finally, for the velocity perturbations, 9 = Vj(5f7*), we 
suppose that they are zero today, since they are decreasing functions and they are precisely zero 
for the strict ACDM model. With this initial ACDM normalization, we may proceed to find the 
spectrum for the RG model in two steps. First, we integrate back the perturbed equations in 
the ACDM case from the present time {z = 0) until some point z = zq, where zq ^ 1. As we 
said, it is not very important the particular value of zq (say 500, 400 etc) provided it lies well 
after the recombination era, i.e. zq < 1100. In this way we determine the values of the three 
fluctuations {Sm,h,v) at z = zq. Second, we can then use these values as initial conditions for 
the RG case. Indeed, since for zq ^ 1 the cosmological constant does not play an important 
role, we can use the initial conditions defined in this manner to solve the perturbed equations 
( 3.17D -( ^.17 ) in our RG model for arbitrary u. In particular, this provides the output for the 
baryonic density spectrum today, 6m{k), for i/ ^ 0. From here the corresponding matter power 
spectrum ( |4.1| ) can be immediately evaluated. Notice that the comparison with the cosmological 
constant case makes sense, since all the difference between the running cosmological constant and 
the particular case of having a strictly constant A (corresponding to the standard ACDM model) 
lies in the evolution of the Universe for a relatively small value oi z. In fact, the contribution of the 
dark energy component begins to become relevant for small redshifts, say for z < 10 or even less, 
and the running is sizeable if |z^| is not smaller than 10~^. To summarize, we use the well known 
semi-analytic expressions for the standard cosmological constant case (which can be reproduced 
from the particular case v = of our general RG framework) to settle the initial conditions for 
our running cosmological constant scenario (z^ 7^ 0). 

The final baryonic spectrum depends essentially on the following three parameters: the relative 
amounts (with respect to the critical density) of baryonic matter, dark matter and dark energy 
today, i.e. O^, 1^2) A/ ^^'^ respectively. The curvature today is given hy = I — — 
where of course = + ^^j^j- Our aim is to compare the theoretical spectrum with the LSS 
data from the 2dF Galaxy Redshift Survey with the error bars evaluated at the la level. To 
start with we consider the standard ACDM model case (therefore we set the relevant RG parameter 
= in our perturbations equations) and produce some numerical examples for the three spatial 
curvatures. The corresponding results for the open (O^ = 0.6), flat (O^ = 0) and closed (0^ = 
—0.5) cases, ah of them with Q^j^^ = 0.21 and n% = 0.04 (as suggested by primordial 
nucleosynthesis), using the transfer function (O), are depicted in Fig. |^a,b,c respectively. The 
presented plots differ essentially by the amount of dark energy and it can be easily seen that the 
best fit corresponds to the flat geometry case (the plot in Fig. |b), with = 0.75. This was the 
expected result and therefore this preliminary exercise serves us as a good normalization of our 
computation. 

In Fig. we fix the values n% = 0.04 and n\ = 0.75 and vary the dark matter content 
in each plot. Specifically, we plot the three cases ^^£)M ~ 0.01, 0.21 and 0.41, representing an 
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Figure 1: Baryonic power spectrum for the ACDM model, for fixed fi^ = 0.04, = 0.21 and for (a) 

51° = 0.15 (open Universe), (b) 51^ = 0.75 (flat Universe) (b) and (c) fi^ = 1.25 (closed Universe), together 
with the LSS data from the 2dfFGRS B. The ordinate axis represents P{k) — \6m{k)\'^ at 2 = as given 



by (4.1) with the transfer function (4.3) and growth function (4.4), while in the abscissa we have the wave 



number k given in h Mpc ^ units. 



open, flat and closed Universe, respectively. We note that the pressureless components of the 
matter content (dark matter and baryons) tend to increase the power spectrum. One can see (for 
the given values of and Q\) that the best fit occurs for the expected DM content, namely 
around = 0.21 (cf. Fig. ^3). This provides further confidence in the normalization of our 

computation with respect to the standard ACDM case. 

4.3 Power spectrum for the running A case 

Now we can start to deal with our main task, that is, to explore the influence of the RG parameter 
v on the processed power spectrum. It is apparent that the larger is the value of ly, the more 
important is the contribution of the function q{z), Eq. (|3.16| ), in the solution of our set of differential 
equations. This function appears as a factor of the Laplacian acting on the perturbed quantities, 
see Eq. ( p.9[ ). Hence, large q{z) entails a more intensive damping of the perturbations contributing 
to depress the power spectrum. This feature can be seen at work in the plots of Fig.^,b,c,d. 
These plots illustrate the power spectrum for running A in the flat case and for non-vanishing and 
positive u. They correspond to = 0.21, n% = 0.04 and = 0.75, for u = 10-^,10^^, 





Figure 2: As in Fig^ (ACDM model) for fixed ^ 0.04, fl^ = 0.75 and for (a) n%j,.j = 0.11 (open 
Universe), (b) ri^JA/ = 0.21 (flat Universe) and (c) il^JA/ — 0-41 (closed Universe), together with the LSS 
data from the 2dfFGRS §. 



10~^ and 10~^ respectively. We can easily grasp the impact of the u parameter, namely we 
see quite patently that the larger is > the more suppressed is the power, mainly at small 
scales (i.e. at large k). For values of > up to = 10~^ the agreement with the LSS data is 
almost perfect; for ly = 10~^ (not shown in the figures) the theoretical curve starts to depart maily 
in the high k region, i.e. at small distances. Above u = 10~^ the spectrum presents a strong 
deviation (depletion of the power) with respect to the observational data, as can be seen in Fig. ^d 
corresponding to = 10~^. This is a quite robust result, and we can state that values of v above 
10~^ are definitely ruled out by the perturbations analysis. 



The interpretation of this result in terms of particle physics can be obtained from Eq. (2.4), 
in which we recall that M stands for the effective mass of the heavy particles contributing to the 
running of the cosmological term. Therefore, in the framework of the given model for the running 
CC - that is, assuming the possibility of an energy exchange between vacuum and matter sector 
via the conservation equation ( |2.2D ~ the exclusion of values u > 10~^ means that we can rule out 
the existence of particles with masses of the Planck mass order, unless their contributions to the 
CC cancel quite accurately due to some strong symmetry (e.g. supersymmetry). 

On the face of the previous result, the next relevant question is: which is then the minimum 
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Figure 3: As in Fig^ but for tlie z^-dependent power spectrum of tlie running cosmological constant 
model, together with the LSS data from the 2dfFGRS[|j. The ordinate axis represents P{k) — \6m{k)\'^ 
where 5m{k) is the solution of ( ^.17l )-( 3.19| ) at z = 0. In all cases we have fixed = 0.04, fi^JA/ = 0-21 



and n\ = 0.75 with the following (positive) values of v: (a) 10"®, (b) lO'^, (c) 10""* and (d) 10" 



upper bound on u compatible with the observational data? It turns out that for values of z/ > 
between 10^^ and 10~^ (and above all between 10~^ and 10~^) it is not possible to state, at this 
stage, if they are definitely ruled out or not. In principle, from the inspection of Fig. ^ we could 
strongly argue in favor of the safer range ^ ^ 10~^. However, we must take into account 
the potential effect of other components. For instance, if we increase the amount of dark matter 
(which is not excluded), the values of > 10~^ may lead to a better fit. On the other hand, 
for negative values of v the bound becomes essentially stronger, i.e. < values produce larger 
deviations beyond a critical range of |z^|. In practice this means that for u outside the interval 
— 10~^ ^ 1/ ^ the spectrum deviates from observations more strongly than for positive values of 
f outside the range ^ ^ 10~®. This is because the parameter v is multiplied, in the Eq. ( p. 17 ), 
by the /c^-proportional terms which stem from the pressure gradient in the perturbed equations. 
When the parameter v becomes negative, the sign of the corresponding coefficient changes. This 
reflects in the evaluated power spectrum in the following way: while for a positive v the spectrum 
is depressed in the small scale regime, for negative v the spectrum is enhanced in that regime, and 
for a value of v well below —10^^ the spectrum becomes rapidly divergent there. This behaviour 
is exemplified in Fig.|^,b,c for the cases u = -10"^, -10-6 and -IQ-"^ respectively. In Fig.| 



we further elaborate on the v < features studied in the previous figure, namely we superimpose 
the detailed evaluation of the spectrum for the five cases u = -10"^, -10"'^, -10"^, -10"^ and 
-10-^. In this way we can better appreciate the aforementioned divergence at small scales, which 
becomes more and more evident the larger is In fact, in this case we see that for \u\ > 10~^ 
the deviation is far more pronounced than for > and this demonstrates that the departure of 
the predicted spectrum from the experimental values in the range -10"^ ^ ^ -10~^ is much 
more violent than in the corresponding positive range 10~^ ^ ^ 10^^. 

The change of behaviour in passing from positive to negative values of ly could be expected. It 
corresponds to the change of an oscillatory behaviour into an exponential one or, in other words, 
to the transition from a harmonic oscillator regime to an anti-harmonic oscillator regime. We can 
also understand why the spectrum is not essentially affected at large scales by the change of sign 
in v, since at large scales the pressure gradient is negligible. For negative the variation of the 
amount of dark matter and dark energy affect the final spectrum in the same way as for a positive 
u case. 

To summarize, the approximate range of values of the RG parameter which are amply allowed 
by the analysis of the matter power spectrum of perturbations (i.e. what we may call the "safest 
range" of allowed values of v) is the following: 

-10"*^ ^ 1/ < 10"^ . (4.9) 

There are nevertheless some ly intervals outside the safest region which can still be considered 
allowed (or at least not ruled out) by the density perturbations analysis. For instance, positive 
values of u in the range 10~^ — 10~^ (mainly those between 10~^ and 10~^) cannot be definitely 
ruled out; we can only say that the corresponding spectrum starts to show some deviation from 
the LSS data becoming more and more acute the closer we get to 10~^. For ly < the exclusion 
effect is stronger, and indeed we can safely rule out the range u < —10^^ for which the predicted 
spectrum becomes rapidly divergent at low scales. Furthermore, for > 10~^ the deviation is 
markedly patent for any sign of ly and we can assert that this range is manifestly excluded. 

Finally, we may ask ourselves how to better improve the aforesaid limits on the RG parameter 
ly. In our opinion this can only be accomplished by the simultaneous crossing of the LSS data with 
further observational data, most significantly with the wealth of CMB anisotropy data. This should 
constraint much better the value of ly beyond the LSS limits that we have been able to obtain. In 
particular, it would probably help to decide what is the situation with the unsettled positive range 
u = 10^^ — 10~^ where the LSS data cannot give a last word. In principle, the CMB spectrum 
predicts an almost spatially flat Universe. But this result is model dependent and another full 
fledged direct evaluation of the cosmological parameters using the running cosmological constant 
model must be made. We plan to perform this analysis in the future. 

5. Discussion and conclusions 



The possibility that the cosmological constant A can be a running parameter within QFT in 
curved space-time was investigated extensively in [ 11 1 . The phenomenological implications of these 
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Figure 4: As in Fig^ (running cosniological constant model), but for fixed fl^ ~ 0.04, fl%j^f — 0.21 and 
— 0.75 and for the following (negative) values of i/: (a) —10^^, (b) —10^®, and (c) —10^^, together with 
the LSS data from the 2dfFGRS il. 



renormalization group (RG) cosmologies on supernovae observations have been amply addressed 
in different works preceding this one p^, IsO, 31, |3^, I 
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Here we have extended the 



analysis of phenomenological implications to the very relevant issue of structure formation. We 
have studied in detail the power spectrum of matter density perturbations for the FLRW type of 
cosmological model with running A. This model is characterized by a fundamental parameter 
associated to the running of A. At the same time this parameter acts, in the present realization of 
the model, as a coupling between the vacuum energy and matter (mainly CDM). There are other 
implementations of this framework where does not lead to this kind of coupling, see |3^. However, 
in the present work we have just focused on the original RG model [ p9| , |3^ , where this coupling is 
present and, therefore, the cosmological term can decay into matter and vice versa. Then, under 
the simplifying assumptions of adiabaticity (i.e. of having no perturbations associated to entropy 
exchange) and of having no anisotropic stress contributions, we have derived the complete coupled 
set of matter density and metric perturbations for this model. The RG running of A is based on 



the assumption of a standard quadratic decoupling law at low energies [29, |3^ and we have taken 
the simplest covariant form of this equation - see (|3.3| ) - as our ansatz to derive the perturbations. 
This phenomenological input leads to a cosmological model with the if^-dependence of the vacuum 
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Figure 5: As in Fig^ (running cosmological constant model), but for fixed = 0.04, ri^JAf = 0-21 a-nd 
= 0.75 with different negative values of v indicated in the figure, together with the LSS data from the 
2diFGRS . The curves corresponding to the last two values of v in the list almost coincide. 



energy density [|Tl| . It is important to emphasize that this does not mean that A ~ tP' , but rather 
that the (renormalization group) variation of A satisfies bh. ^ H^, and hence A = a + bH^ - see 
Eq. (^.5|). The effect produced by the running depends on a single parameter u, which characterizes 
the particle mass spectrum of the quantum theory below the Planck scale. For = we retrieve 
the standard ACDM model with time-independent cosmological constant, while the values around 
the "canonical" value v = ^ 10~^ correspond to having effective particle contributions of order 
of the Planck mass. 

The numerical analysis of the perturbations proves to be in qualitative agreement with the 
results of the previous estimate in Ref. |^^, based on restraining the amplification of the density 
fluctuating spectrum at the recombination era caused by the vacuum decay into CDM. In both 
cases (viz. the direct calculation and the previous estimate) it is found that the "canonical" 
value of u leads to strong deviations from the observational data. For positive values of u in 



X 



the range 10~° ^ ^ 10~ , corresponding approximately to the GUT particle spectrum (M 
(10~^ — 10^'^) Mp), we meet relatively small deviations of the theoretical curve from the = 
case. The deviations are higher the larger is |z^|. Let us recall that the upper bound on v obtained 
from the phenomenological estimate of [40| is 10^'^ and applies only for u ^ 0. This bound is not 
very far away from our own most conservative result \u\ < 10~^, although in our case it applies 
for both signs of u. Moreover, the approximate bound is much weaker than the one defined by the 
so-called "safest range" - cf. Eq. ([4.9[). However, as we already emphasized, we cannot definitely 



exclude the u > interval 10 ^ — 10 ^ (and specially the subinterval 10 ^ — 10 ^). Therefore, we 
may assert, in the most possible conservative way, that our direct bounds on u > are between 



one to two orders of magnitude more restrictive than those of [40|, but of course they have been 
derived in a much more rigorous way. Let us also remark that our exclusion of the possibility that 
u could be of the order of magnitude of the "canonical" value i^o ~ 10^^ is at variance (by at 
least two orders of magnitude) with the estimates obtained from the approximate methods used 
in references |41, 42| 



It is also interesting to compare our results with those of Ref. [47|, where a detailed study 
of matter density perturbations in interacting quintessence models is addressed. These authors 
perform their study also under the assumption of an adiabatic regime and with no contributions 



from anisotropic stress. In this kind of modified quintessence models, first proposed in |54], 
one couples matter to quintessence by means of a source function Q which depends on a free 
parameter, called cP in Ref. ||4^. This is the only free parameter in that model, similarly as 
is I' in our own framework under the specified common set of assumptions. In the models of 
interactive quintessence, the dark energy (DE) can also decay into matter (parallel to our case 
where A can also decay into matter) and as a result the amount of DE in the past is larger than 
in ordinary quintessence models with self-conserved DE. This corresponds, in our framework, to 
have V > because then A becomes larger in the past. The authors of [^] find that the effect of 
the c^-coupling is to introduce a damping of the power spectrum of matter perturbations at low 
scales; that is to say, the growth of density perturbations for large wave numbers becomes smaller 
as compared to quintessence models with uncoupled DE to matter. In our RG framework, we 
observed a similar depletion of the spectrum at low scales for v > as compared to = (the 



standard ACDM model) - cf. sections O and O. Moreover, we have found numerical bounds 
on V that are quantitatively similar, in fact a bit more stringent than those on obtained in the 



model of [47|. These authors found that their power spectrum is consistent with the 2dFGRS data 
provided < 10~^, and we find that the largest allowed value of u by the same LSS data is of 
order 10~^. The basic agreement among the different analyses suggest that the loss of power at 
low scales could be a general feature of dynamical models in which the dark energy decays into 



matter during the cosmic evolution. The latter situation is always the case in Ref. |47] because 
> 0, and it is also our case for > 0. However, in our RG framework, z/ can have any sign 
(depending on the dominance of bosons or fermions in the loop contributions to the running of 
A) and so we have an opportunity to explore also the case z/ < 0. For this sign of u, the running 
A becomes progressively smaller in the past and reaches eventually large negative values (see the 
detailed evolution plots of Ref. [^). As a result we may expect an overproduction of structure 
associated to the fact that a large negative A increases the gravitational collapse. Since, however. 



this is not observed (cf. the plots of section 4.3 detecting an excess of power at low scales for 



z^ < 0), this may provide a physical explanation of why negative values of z^ are significantly more 
restrained. While the safest and most stringent range of values of is expressed in a nutshell 
in Eq. ([4.91), we cannot exclude that could be higher (mainly in the z^ > case) once it will 



^The results of these references were presented in terms of a parameter e, and they concluded that values of e at 
the level of 10"^ are allowed. However, e is nothing but 3 times the parameter u originally defined in |29]. 



be possible to combine the analysis of the matter power spectrum with the corresponding CMB 
analysis. At the same time, we should recall that the bounds obtained on v in our case, and on (? 



in 1 47 1, could suffer some renormalization if we would abandon the set of simplifying hypotheses 
that have been made. However, our study shows at least the kind of impact to be expected on 
the matter power spectrum predicted by cosmological models with running A in interaction with 
matter. Further investigations would be needed to assess the effect in a more thorough way, but 
they lie beyond the scope of the present work. 

Let us also finally remark that the bounds we have placed on the RG parameter within our 
set of assumptions, can be interpreted on purely phenomenological grounds, namely as bounds on 
a parameter that gauges the possible degree of interaction between CDM and the vacuum energy, 
even if no fundamental RG model is invoked. In principle, the phenomenological implications of 
a non- vanishing v could be potentially detectable in the future supernovae observations, either 
through a direct measure of the running of A []32| , |33[| or (more efficiently in practice) through the 
possibility to observe the redshift evolution of the effective equation of state parameter associated 



to the variable CC model [35, 36|. However, the very study of the evolution of the density and 
metric perturbations, which has led us to firmly exclude the range \v\> 10~^, already demonstrates 
the difficulty of the previous methods and at the same time emphasizes the great advantage of 
the density perturbations approach. Indeed, the strong sensitivity of the matter power spectrum 
to a variable cosmological term interacting with the CDM shows that this is perhaps the ideal 
observable to look at. Therefore, future data on LSS could eventually provide the clue to unravel 
an interaction (even if very small) between matter (essentially CDM) and vacuum, and at the 
same time it may hint at the possible RG physics behind it. Let us conclude by mentioning that it 
should be very interesting to expand our analysis to the aforementioned second model of running 



CC [34|. In this alternate model the cosmological term is also evolving in time and redshift, but 
there is no energy exchange between matter and vacuum and we may expect much more freedom 
in choosing the value of v. We leave this study for a future publication. 
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